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Abstract 

In this paper we study the spectrum of certain large random Hermitian Jacobi matrices. These 
matrices are known to describe certain communication setups. In particular we are interested in an 
uplink cellular channel which models mobile users experiencing a soft-handoff situation under joint 
multicell decoding. Considering rather general fading statistics we provide a closed form expression 
for the per-cell sum-rate of this channel in high-SNR, when an intra-cell TDMA protocol is employed. 
Since the matrices of interest are tridiagonal, their eigenvectors can be considered as sequences with 
second order linear recurrence. Therefore, the problem is reduced to the study of the exponential growth 
of products of two by two matrices. For the case where K users are simultaneously active in each cell, 
we obtain a series of lower and upper bound on the high-SNR power offset of the per-cell sum-rate, 
which are considerably tighter than previously known bounds. 
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I. Introduction 

The growing demand for ubiquitous access to high-data rate services, has produced a huge 
amount of research analyzing the performance of wireless communications systems. Cellular 
systems are of major interest as the most common method for providing continuous services to 
mobile users, in both indoor and outdoor environments. Techniques for providing better service 
and coverage in cellular mobile communications are currently being investigated by industry and 
academia. In particular, the use of joint multi-cell processing (MCP), which allows the base- 
stations (BSs) to jointly process their signals, equivalently creating a distributed antenna array, 
has been identified as a key tool for enhancing system performance (see [1][2] and references 
therein for surveys of recent results on multi-cell processing). 

Most of the works on the uplink channel of cellular systems deal with a single-cell setup. 
References that consider multi-cell scenarios tend to adopt complex multi-cell system models 
which render analytical treatment extremely hard (if not, impossible). Indeed, most of the results 
reported in these works are derived via intensive numerical calculations which provide little 
insight into the behavior of the system performance as a function of various key parameters 
(e.g. [3]-[8]). 

Motivated by the fact that mobiles users in a cellular system "see" only a small number of 
BSs, and by the desire to provide analytical results, an attractive analytically tractable model for 
a multi-cell system was suggested by Wyner in [9] (see also [4] for an earlier relevant work). In 
this model, the system's cells are ordered in either an infinite linear array, or in the familiar two- 
dimensional hexagonal pattern (also infinite). It is assumed that only adjacent-cell interference 
is present and characterized by a single parameter, a scaling factor a e [0, 1]. Considering non- 
fading channels and a "wideband" (WB) transmission scheme, where all bandwidth is available 
for coding (as opposed to random spreading), the throughput obtained with optimum and linear 
MMSE joint processing of the received signals from all cell-sites are derived. Since it was first 
presented in [9], "Wyner- like" models have provided a framework for many works analyzing 
various transmission schemes in both the uplink and downlink channels (see [2] and references 
therein). 

In this work we consider a simple "Wyner-like" cellular setup presented in [10] (see also 
[11]). According to this setup, the cells are arranged on a circle (or a line), and the mobile users 
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"see" only the two BSs which are located on their cell's boundaries. All the BSs are assumed 
to be connected through an ideal backhaul network to a central multi-cell processor (MCP), 
that can jointly process the uplink received signals of all cell-sites, as well as pre-process the 
signals to be transmitted by all cell-sites in the downlink channel. The users are hence in what 
is referred to as a "soft-handoff ' situation, which is very common in practical real-life cellular 
systems, and is therefore of real practical as well as theoretical interest (see for example [12] 
for a recent survey on handoff schemes). With simplicity and analytical tractability in mind, 
and in a similar manner to previous work, the model provides perhaps the simplest framework 
for a soft-handoff setting in a cellular system, that still represents real-life phenomena such as 
intercell interference and fading. 

Unfortunately, the analysis of "Wyner-like" models in general and the "soft-handoff setup in 
particular presents some analytical difficulties (see Section III-B|) when fading is present. These 
difficulties render conventional analysis methods such as large random matrix theory impractical. 
Indeed the per-cell sum-rate rates supported by MCP in the uplink channel of the "soft-handoff 
setups are known only for limited scenarios such as non-fading channels, phase-fading channels, 
fading channels but with large number of users per-cell, and Rayleigh fading channels with single 
user active per-cell [10] [1 1] [13] . The latter result is due to a remarkable early work by Narula 
[14] dealing with the capacity of a two-tap time variant ISI channel. Calculating the per-cell 
sum-rate capacity supported by the uplink channel of the "soft-handoff setup in the presence 
of general fading channels (not necessarily Rayleigh fading channels), when finite number of 
users are active simultaneously in each cell remains an open problem (see [1 1 ] [1 5] for bounds 
on this rate). As will be shown in the sequel, this problem is closely related to calculating the 
spectrum of certain large random Hermitian Jacobi matrices. The high-SNR characterization of 
the sum-rate capacity, previously unknown, is the main focus of this work. 

In particular we calculate the high-SNR slope and power offset of the rate with a single user 
active per-cell (intra-cell TDMA) under a rather generic fading distribution. We also prove the 
following results for any given number of active users per-cell. We prove the existence of a 
limiting sum-rate capacity when the number of cells goes to infinity and calculate the high-SNR 
slope in Theorem [2l Moreover, we give bounds on the high-SNR power offset in Proposition 
[U In particular, we give a sequence of explicit upper- and lower-bounds; the gap between the 
lower and the upper bounds is decreasing with the bounds' order and complexity. 
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The rest of the paper is organized as follows. In Section [n] we present the problem statement 
and main results. Section [III] includes a comprehensive review of previous works. Several ap- 
plications of the main result are discussed in Section [IV] Concluding remarks are included in 
Section [V] Various derivations and proofs are deferred to the Appendices. 

II. Problem Statement and Main Results 

A. System Model 

In this paper we consider a linear version of the cellular "soft-handoff ' setup introduced in 
[10] [11], according to which M + 1 cells with K single antenna users per cell are arranged on 
a line, where the M single antenna BSs are located on the boundaries of the cells (see Fig. Q] 
for the special case of M = 3). Starting with the WB transmission scheme where all bandwidth 
is devoted for coding and all K users are transmitting simultaneously each with average power 
p, and assuming synchronized communication, a vector baseband representation of the signals 
received at the system's BSs is given for an arbitrary time index by 

y = H M x + n (1) 

The M x K(M + 1) channel transfer matrix Hm is a two block diagonal matrix defined by 

( a x &i ••■ ^ 



M 







(2) 



: 
\ • ■ • a M b M J 

where a m and b m are 1 x K row vectors denoting the channel complex fading coefficients, 
experienced by the K users of the mth and (m + l)th cells, respectively, when received by the 
mth BS antenna, n represents the M x M zero mean circularly symmetric Gaussian noise vector 
n~CM{Q,I M ). 

We assume throughout that the fading processes are i.i.d. among different users and BSs, with 
a>m,k ~ an d b m k ~ 7r b , and can be viewed for each user as ergodic processes with respect 
to the time index. We denote by P the probability associated with those random sequences and 
by E the associated expectation. We will be working throughout with a subset of the following 
assumptions. 
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(HI) E^ a (log \x\f < ocQ and E nb (log \x\) 2 < oo. 

(H2) 7r a and 7r& are absolutely continuous with respect to Lebesgue measure on C. 

(H3) There exists a real M. such that if x is distributed according to n a (resp. n^) then the 

density of |x| 2 is strictly positive on the interval [Ai; oo). 
(H3') There exist m a < M a G R + U {00} (resp. m 6 < .M& G M + U {00}) such that if x is 

distributed according to n a (resp. n^) then the density of \x\ and the Lebesgue-measure 

on [m a ; Ai a ] (resp. [m fe ; A^J) are mutually absolutely continuous. 
(H4) There exists a ball in C such that the Lebesgue measure outside that ball is absolutely 

continuous with respect to % a and Tif,. 
We further assume that the channel state information (CSI) is available to the MCP only, 
while the transmitters know only the channel statistics, and cannot cooperate their transmissions 
in any way. Therefore, independent zero mean circularly symmetric Gaussian codebooks conform 



with the capacity achieving statistics, where x denotes the (M + 1)K x 1 transmit vector x ~ 
CA/"(0, pImk), and p is the average transmit power of each user|^ (p is thus equal to the transmit 
SNR of the users). 

With the above assumptions, the system © is a multiple access channel (MAC). We are 
interested in the per-cell sum-rate capacity 

Cm (P) — (log det Gm) [nats/channel use] , (3) 

where P = Kp is the per-cell transmitted average power, 

Gm — I m + pH h[H M , (4) 

and the expectation is taken over the channel transfer matrix entries. (Here and in the sequel, for 
a scalar z G C, denotes the complex conjugate, while for a matrix A, denotes the matrix 

'A natural base logarithm is used throughout this work unless explicitly denoted otherwise. 

2 Note that since the channel transfer matrix Hm is a column-regular gain matrix (see definition in [16]) when M — * oo, the 
capacity achieving statistics remains the same in this case, even if we allow the users to cooperate as long as they are unaware 
of the CSI. 
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with A^(i,j) = A(j,iy.) The non-zero entries of the Hermitian Jacobi matrix Gm are equal to 

l G ^]m,m~l = P < b m-l, a m > , 

[ G M\ m , m = 1 + P (l«m| 2 + \b m \ 2 ) , ( 5 ) 
[ G M]m,m+l = P < a m+l] b m > , 

where out-of-range indices should be ignored, and for any two arbitrary L length vectors a, b 

we define < a;b >= Yla=i anc ^ \ a \ 2 ~ < a ' a > - 

Since we shall focus on the asymptotes of infinite number of cells M — > oo, boundary effects 
can be neglected and symmetry implies that the rate © equals the maximum equal rate (or 
symmetric capacity) supported by the channel [17]. 

The above description relates to the WB protocol where all users transmit simultaneously. 
According to the intra-cell TDMA protocol only one user is simultaneously active per-cell, 
transmitting 1/K of the time using the total cell transmit power P. In this case it is easily 
verified that with no loss of generality, we can consider a single user per cell in terms of the 
per-cell sum-rate, setting K = 1 in (OQ) and ©. 

B. Analysis Difficulty 

Many recent studies have analyzed the rates of various channels using results from (large) 
random matrix theory (see [18] for a recent review). In those cases, the number of random 
variables involved is of the order of the number of elements in the matrix Gm (or Hm), and 
self-averaging is strong enough to ensure convergence of the empirical measure of eigenvalues, 
and to derive equations for the limit (or its Stieltjes transform). In particular, this is the case if 
the normalized continuous power profile of Hm, which is defined as 

P „M)4 E( |[H„] ( /) ; ±< r <^, w ± M <t< T l±L R , ,6) 

converges uniformly to a bounded, piecewise continuous function as M — > oo, see e.g. [18, 
Theorem 2.50] and [19] for fluctuation results. In the case under consideration here, it is easy to 
verify that for K fixed, Vm{?, t) does not converge uniformly, and other techniques are required. 

C. Extreme SNR Regime Characterization 

As mentioned earlier, the per-cell sum-rate capacity of the "soft-handoff ' setup is known only 
for certain limited cases to be elaborated in the next section, and in general analytical results are 
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hard to derive. As an alternative to deriving exact analytical results we focus here on extracting 
parameters which characterize the channel rate under extreme SNR scenarios. The reader is 
referred to [20] - [22] for an elaboration on the extreme SNR characterization. 

a) The Low-SNR Regime: This regime is usually the operating regime for wide-band 
systems [21]. 

The average per-cell spectral efficiency in bits/sec/Hz, expressed as a function of the sys- 
tem average transmit SNR, E b /N , is evaluated by solving the implicit equation obtained by 
substituting 



(V) 

N J N 

in ©, where CM(Eb/No) = Cm(-P)/ log 2 stands for the uplink spectral efficiency measured in 
[bits/sec/Hz]. The low-SNR regime is characterized through the minimum transmit E b /N that 
enables reliable communications, 



N 0a 

and the low-SNR spectral efficiency slope 



log 2 

CW(0) 



(8) 



C M (0) 

So " 

-Cm(0) 

yielding the following low-SNR affine approximation 

E b 



C 



M 



No 





(E b 




E b 


/ 0|dB 


\N 


dB 


Nn ■ 

1 ' U mm 



dB 



[bits/sec/Hz]. 



(9) 



(10) 



In the above definitions 3|dB = 101og 10 2, and C*m(0) and Cm(0) are the first and second 
derivatives (whenever exist) with respect to P of the per-cell sum-rate capacity, respectively, 
evaluated at P = 0. Focusing on Gaussian channels with receiver CSI only, it can be shown [21] 
that there is no need to calculate the two derivatives of the rate in P = 0, and that the low-SNR 
parameters are simply given by 



El 
n 0i 



MK log 2 
U:[Erf M H M 



So 



tr ( EH j M H M 



M 



t r I 2 ( H ] M H M 



(11) 
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b) The High-SNR Regime: This is usually the operating regime for high-data rate (high 
spectral efficiency) systems (that is the case actually in all 2.5/3 G standards). 

The high-SNR regime is characterized through the high-SNR slope (also referred to as the 
"multiplexing gain", or "pre-log") 

Soo 4 lim = lim PC M (P) , (12) 

P->oo logP P^oo 

and the high-SNR power offset 

£«, 4 lim -L (logP - ^^y\ , (13) 
yielding the following affine capacity approximation 

c u (P) « %^ (PU - 3| dB £oc) • (14) 

Note that the high-SNR approximation reference channel here is that of a single isolated cell, 
with no fading, and total average transmit power P. 

The high-SNR characterization of the per-cell sum-rate supported by the "soft-handoff ' uplink 
channel is known only in certain limited scenarios (see Section [Til]) and is the main focus of 
this work. 

D. Main Results 

Recall the definition of Cm(P), c.f. ©. Starting with intra-cell TDMA scheme where only 
one user is active per-cell transmitting with power P we have the following. 

Theorem 1 [intra-cell TDMA scheme K — 1, high-SNR characterization] Assume (HI) and 
(H2). 

a) For every P > 0, Cm(P) converges as M goes to infinity. We call the limit C(P). 

b) We get the following bounds on C (P), 

max(E 7ra log(l + P \x\ 2 ),E nb log(l + P \y\ 2 )) < C(P) < E w log(l + P(|a;| 2 + \y\ 2 )). 

c) Further assume [(H3) or (H3')]. As P goes to infinity, 

C(P) = logP + 2max(E Wa log \x\ , log |ar|) + o(l). 
In particular, = 1 and = —2 max (E^ log 2 |x| , E^ log 2 |x|). 
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Note that point c) shows that the lower bound of point b) is tight in the high-SNR regime. 
Proof: The proof of points a) and c) follows from Theorem \5\ of Appendix [A] where we prove 
that the variable Cm(P) — 1/Mlogdet converges almost surely. Note however that 

1 1 - 

< — logdet G M < — lo § (! + P(\ a m\ 2 + K\ 2 )) , (15) 

m=l 

and the second inequality is due to Hadamard's inequality for semi-positive definite (SPD) 
hermitian matrices. With (HI), it follows that Cm(P) is uniformly integrable, and hence the 
almost sure convergence implies convergence in expectation. Recalling that C M {P) = E C M (P) 
completes the proof of point a) and c). 

Let us show point b) using the tools of [23]. We first show the lower bound. We consider n, 
x and y as in ©. 

C M (P) = j-fl (sb; y\{ai)i<i<M, {k)i<i<M) 



1 M 

— ^ I [Xj] y\(Xi)i<i<j, (Oi)i< 4 <M, (bi)i<i< M ) 



3=1 

M 



— ^ %-ll( a '«)l<i<J' ( a i)l<i<M, (i>i)l<i<Af) 

3=1 

1 M 

which is the per-cell sum-rate capacity of a single user fading channel. Therefore, the lower 
bound is [24] log(l + P\y\ 2 ). As argued in the proof of Theorem [5] in Appendix El we 
can exchange the role of ir a and 7r&, thereby getting the claimed lower bound. Finally, the upper 
bound of b) follows immediately from Hadamard's inequality for SPD hermitian matrices. ■ 
In the proof of Theorem [5] (intra-cell TDMA scheme), we use ideas from the theory of 
product of random matrices. Note that Cm{P) = Sm=i 1°§(1 + P\n) where {A m }^f =1 
are the eigenvalues of HmH\ { , and the analysis of capacity hinges upon the study of spectral 
properties of H mH^ m . The main idea is to link the spectral properties of the latter matrix 
with the exponential growth of the elements of its eigenvectors. Since HmH' m is a Hermitian 
Jacobi matrix, hence tridiagonal, its eigenvectors can be considered as sequences with second 
order linear recurrence. Therefore, the problem boils down to the study of the exponential growth 
of products of two by two matrices. This is closely related to the evaluation of the top Lyapunov 
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exponent of the product; The explicit link between Cm(P) an d the top Lyapunov exponent is the 
Thouless formula (see [25] or [26]), a version of which we prove in Appendix iDl We emphasize 
however that we do not use the Thouless formula or Lyapunov exponents explicitly in the proof 
of Theorem |5] 

Like in the result of Narula [14] described below in SectionUm our approach uses the analysis 
of a certain Markov Chain. Unlike [14], we are not able to explicitly evaluate the invariant 
measure of this chain. Instead, we use the theory of Harris chains to both prove convergence 
and continuity results for the chain. The appropriate definitions are introduced in the course of 
proving Theorem \5\ 

We remark that Theorem Q] continues to hold in a real setup, that is if instead of (H2), we 
assume 

(H2') ir a and n b are supported on R and are absolutely continuous with respect to Lebesgue 
measure on R. 

Since the argument is identical, we do not discuss this case further. It is also noted that unlike 
the non-fading case, where intra-cell TDMA scheme is optimal (see [9]), it is proved to be 
suboptimal for K > 1 in the presence of fading [27], yet TDMA it is one of the most common 
access protocols in cellular systems. 

Turning to the WB scheme (which is the capacity achieving scheme [27]), where all the 
bandwidth is used for coding, and all K users are transmitting simultaneously with average 
power p (and total cell average power P = Kp), we have the following less explicit high-SNR 
characterization. 

Theorem 2 [WB scheme K > 1, high-SNR characterization] Assume (HI), (H2) and (H4), and 
K > 1. 

a) For every P > 0, Cm{P) converges as M goes to infinity. We call the limit C(P). 

b) We get the following bounds on C(P), 

max(Elog(l + P|a| 2 /^),Elog(l + P|b| 2 /^)) < C{P) < Elog(l + P(\a\ 2 + \b\ 2 )/K), 

where the expectation is taken in the following way: the random variables a and b are 
independent, and a (resp. b) is a complex K -vector whose coefficients are independent 
and distributed according to n a (resp. n b ). 
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c) As P goes to infinity, 



C(P)=\ogP + E\og\ 6 -^- ) +o(l) , (16) 



where the expectation is taken in the following way: the random variables e and b 
are independent, and b is a complex K-vector whose coefficients are independent and 
distributed according to 7r b . The law of e is m , which is the unique invariant probability 
of the Markov chain defined by 

2 / e n + |b n _i| 2 sin 2 (a n ,6 n _i) \ 

e n+1 = \a n \ -2 , (17) 

V e n + |6 n _i| J 

where for any two arbitrary equal length vectors a, b, 

sm 2 (a,b)±l- l<a ^ b> f . (18) 
I a. I 161 



In particular, 5^ = 1 and = — Elog 2 ( e+ ^ \ 



As with the case K — 1, point a) and c) of Theorem [2] follow from the almost sure convergence 
stated in Theorem |2T1 of Appendix O using (HI) and (fl"5l) . As with Theorem |5l we do not use 
the Thouless formula or Lyapunov exponents explicitly in the proof of Theorem |2U The proof 
of point b) is the same as the proof of Theorem [TJb). It is worth mentioning that in contrast to 
Theorem [H the non-asymptotic lower bound b) is not tight in general for large SNR. This is 
since it is an increasing function of K and converges to a rate of a single-user Gaussian scalar 
channel, which is smaller than the asymptotic rate of (|23|) . 

Note that although the roles of the sequences {a n } and {6„} in (fTTT) are not symmetric, the 
expression (fT6l) is symmetric in n a and 7Tb, as is the case for K = 1. 

We conclude this section by noting that while Theorem [2] (WB scheme K > 1) does not give 
explicit expressions for the high-SNR power offset as Theorem [Q its proof leads immediately to 
easily computable bounds. In the following, the notation is as in Theorem |2l and we let e n (a) 
denote the Markov chain (fT71) . with initial condition e (a) = a. 

Proposition 3 Assume (HI), (H2) and (H4), and K > 1. Then, 

Eiogi < |c( p) _ logP] < Elog to + i"i 2 ' 
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where the expectation is taken in the following way. e„(0) (resp. e n (oo)) and b are independent, 
b is a complex K-vector whose coefficients are independent and distributed according to 7r&. 
e n (0) (resp. e n (oo)) is the n-th step of the Markov chain defined by rfi"7l) with initial condition 
e (0) = (resp. e (oo) = ooj. 

Indeed, since the expression (fTTT) for e n+ i is monotone increasing in e n , the law of e in Theorem 
[2] is stochastically dominated below by the law of e n with intial condition 0, and stochastically 
dominated above by the law of e n with initial condition oo. That same monotonicity also shows 
that the sequences of laws of e n (0) (resp., e n (oo)) are monotone increasing (resp., decreasing) 
with respect to stochastic order. 

As a direct consequence of Proposition [3] with n = 1 and (TT3l) . we get the following bounds 
on the high-SNR power offset 

_ f \a\ 2 S m 2 (a,b) + \b'\ 2 \ 
R . , < Coo < -E log 2 I i-! K J > ( 19 ) 

where the expectation is taken in the following way: a, b and b' are independent, and a (resp. 
b, b') is a complex if- vector whose coefficients are independent and distributed according to n a 
(resp. 7Tb). Note that for K going to infinity, if we assume ix a = 7r& and zero mean, then sin 2 (a, b) 
converges to 1, therefore the ratio between the upper- and lower-bound of (19), converges to 1, 
which also agrees with the asymptotic result of (1371) . 

Numerical Results: In Figures [2] and [3] we present the high-SNR power offset bounds of 
Proposition [3] in the special case of Rayleigh fading (real and imaginary parts are independent 
Gaussian random variables with zero mean and variance l/y/2), for K = 2 and K = 10 users 
per-cell respectively. The curves are produced by Monte Carlo simulation with 10 5 samples. The 
figures include also the lower bound of [11], see (|38l) . and the asymptotic results (and lower 
bound) for large number of users per-cell = — 1 (achieved by taking K to infinity in (1381)). 
Examining the figures it is observed that the new bounds are getting tighter with their order n 
and that the new lower bound is tighter than (1381) already for n = 2. Moreover, fixing the order 
n, the new bounds are getting tighter with the number of users per-cell K. This observation 
is also evident from Fig. 51 where the bounds are plotted for a fixed order n = 2 versus the 
number of users per-cell K. Finally, since the upper bound of Fig. [2] is negative, we conclude 
that the presence of Rayleigh fading is beneficial over non-fading channels in the high-SNR 
region already for K = 2. (See [11] for a similar conclusion in the low-SNR region.) 
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III. Background, Previous Results and Bounds 

In this section we briefly summarize previous work on the "soft-handoff ' uplink cellular model 
introduced in [10] [1 1] . For conciseness, we restrict the discussion to the case where ix a = 7r&. 
Most of the results in the sequel can be extended to include the general case where ir a ^ n h . 

Starting with non-fading channels (i.e., when 7r a and ix^ are singletons at 1), the per-cell 
sum-rate capacity of the uplink channel is given for M — > oo by [11] 



Rat = log I I . (20) 

This rate is achieved by any symmetric intra-cell protocol with average transmit power of P 
(e.g. intra-cell TDMA, and WB protocols). It is noted that the same result holds also for phase 
fading processes [13]. 

The extreme SNR characterization of (|20l) is summarized for the non-fading setup by 

So = - , -rf- = f , 5*00 = 1, Coo = . (21) 

■J ^ V min * 

Returning to the flat fading setup, the channel coefficients are taken as i.i.d. random variables, 
denoting by 

m 1 = E(a m>k ) = E(b m>k ) ; m 2 = E(|a mifc | 2 ) = E(\b rn:k \ 2 ) 

m A > V m,k (22) 

m 4 = E(|a mifc | 4 ) = E(\b m , k f) ; K = — ^ 

the mean, second power moment, fourth power moment and the kurtosis of an individual fading 
coefficient. 

The per-cell sum-rate capacity of the WB scheme with fixed P and increasing number of 
users and cells M, K — ■> oo, is given by [llj^l 



1 + 2Pm 2 + \ 1 + APm 2 + 4P 2 (m 2 2 - 
iCb-f = log I y - : — - I . (23) 




The rate is maximized for a zero mean fading distribution and is given by 

i? wb _ f = log(l + 2m 2 P) . (24) 

3 Here, the number of users K is taken to infinity and then the number of cells M is taken to infinity. 
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Comparing (|20l) and (|24l) (with m 2 = 1), it follows that the presence of fading is beneficial in 
case the number of users is large. We note that (|23l) is also shown in [11] to upper bound the 
respective rate for any finite number of users K. 

Returning to the intra-cell TDMA (K = 1), for which standard random matrix theory is not 
suitable (see Sec. III-BI) . the powerful moment bounding technique employed in [27] for the 
Wyner model, can be utilized to obtain lower and upper bounds on the per-cell sum-rate. 

An alternative approach which replaces the role of the singular values with the diagonal 
elements of the Cholesky decomposition of the the matrix Gm, was presented by Narula [14] 
for a two diagonal nonzero channel matrix Hm whose entries are i.i.d. zero-mean complex 
Gaussian (Rayleigh fading). Originally, Narula had studied the capacity of a time varying two 
taps inter-symbol-interference (ISI) channel, where the channel coefficients are i.i.d. zero-mean 
complex Gaussian. With the above assumptions regarding the ISI channel coefficients it is easy 
to verify that the capacity of this model is equal to the per-cell sum-rate capacity of an uplink 
intra-cell TDMA scheme employed in the "soft-handoff" model. 

Following [14], we use the Cholesky decomposition applied to the co variance matrix of the 
uplink intra-cell TDMA scheme output vector Gm = L m DmUm, where L M (resp. U m) is a 
lower triangular (resp. upper triangular) matrix with 1 on the diagonal. The diagonal entries of 
Gm are given (with K = 1) by 

d m = l + P\a m \ 2 + P\b m \ 2 (l-P^p^-) , m = 2,...,M, (25) 

y / 

where the initial condition of (1251) is d\ = 1 +P \a\\ 2 + P \bi\ 2 . Thus, the diagonal entries {d m } 
form a discrete-time continuous space Markov chain; Narula's main observation was that this 
chain possesses a unique ergodic stationary distribution, given by 

log(ar)e-f 
hl Kp) f 

where Ei(x) = exp [~^ dt is the exponential integral function. Further, as is proved in [14], 
the strong law of large numbers (SLLN) holds for the sequence {log<i m } as M — > oo. Hence, 
the average per-cell sum-rate capacity of the intra-cell TDMA scheme (K = 1) can be expressed 
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as 

■Rtdma-f = lim E ( 77 lo S det G M 



Inn 2 ( — log det (L M D M U M ) J (27) 



M^oo 



/ 1 M \ 

\ m=0 / 



M- 

where the last expectation is taken with respect to fd(x), as defined in (|26l) . In particular, 

/« (log(x)) 2 e-f 

-Ktdma-f — / /n . {16) 

Narula's approach is based on an explicit calculation of the invariant distribution f d , and is thus 
tied to Rayleigh fading. Modifications of key parameters (such as the entries' PDF, and the 
number of nonzero diagonals) lead to analytically intractable expressions. 

Another result derived by following the footsteps of [14] is an upper bound on the per-cell 
sum-rate of the WB scheme with finite K and infinite number of cells M — > oo, in the presence 
of a general fading distribution, given by 



1 + 2Pm 2 + Jl + APm 2 + 4P 2 (l - [m\ - |mj| ) 
Pwbk-f < log I J . (29) 

and in the special case of zero mean unit power (m x = 0, m 2 = 1) fading distribution (e.g. 
Rayleigh fading) the bound reduces to 



Pwbk-f < log I I . (30) 

This result which is proved in [14] for K = 1 (intra-cell TDMA protocol) and expanded to an 
arbitrary K in [15], is derived by noting that the average of the determinant of the received 
vector covariance matrix Gm can be recursively expressed by 

E(det G m ) = A E(det G m _i) - B E(det G m _ 2 ) ; m = 3, . . . , M , (31) 

with initial conditions 

E(det Gi) = A ; E(det G 2 ) = A 2 - B , (32) 

where 

p2 

A=l + 2Pm 2 ; B = — [m\ + (K - 1) Kl 4 ) . (33) 
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See Appendix |E] for more details. The solution to (I3TT) is given by 

E(det G m ) = (f r m - <j) s m , (34) 

where 

r = ~ (A + vM 2 - 45) ; s = ~(A- VA 2 - 45) , (35) 

are real and positive, and cp, <p are determined by the initial conditions (l32l) . Finally, (T29l) is 
derived by the following set of inequalities 

#wbk-f = hm -^-E(logdetG M ) < lim -!- logE (det G M ) = logr , (36) 

M—>oo JVl M^oo iVl 

where the inequality is due to Jensen's inequality, and the last equality follows from the fact that 
r > s, and M — > 00. In the case of K = 1, the upper bound of ([301 ) coincides with the per-cell 
sum-rate capacity of the non-fading setup ([20l . Thus, the presence of Rayleigh fading decreases 
the rates of the intra-cell TDMA protocol supported by the "soft-handoff ' model. Nevertheless, 
it is shown in [11] that already for K = 2 the presence of fading may be beneficial at least 
for low SNR values. The tightness of the bound is demonstrated by noting the for K — > 00 it 
coincides with the asymptotic expression of (|23l) . 

The extreme SNR characterization of the WB rate for M — > 00 in the presence of a general 
fading distribution is summarized by [11] 

S = 2 ■ — = log2 

g 1 l m il 4 ~^]" ' No min 2m 2 

- h 2m 2 . (37) 



log 2 ^m 2 + \J(l-i) (ml - |mx| 4 )^ < £ c 



5*00 ^ 1 ) 

The bounds of the high-SNR parameters are tight for K ^> 1. For the special case of Rayleigh 
fading the extreme SNR characterization are given by [11] 

c __ 2 . Eb _^ log_2 

i+5 ' Nomin 2 



/ / T\ (38) 

Soo = 1 ; -to & (l + ^/l - j?) < £00 < ^ . 

where 7 w 0.5772 is the Euler-Mascheroni constant. It is noted that the right inequality of the 
high-SNR power offset is tight for K = 1, while the left inequality is tight for K ^> 1. The 
beneficial effects of Rayleigh fading and increasing number of users are evident when compared 
to the non-fading extreme-SNR parameters of the respective non-fading setup (I2TI) . 

To conclude this section we emphasize that calculating exact expressions for the high-SNR 
parameters of the WB protocol rate with finite number of users per-cell and general fading 
distribution remains an open problem. 
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IV. Applications 

In this section we present several applications of the main results presented in this work (see 
Section HlDl). 

c) Intra-Cell TDMA and Rayleigh Fading: Assuming that only one user is active per-cell 
K = 1 and symmetric Rayleigh fading channels (i.e. 71^2 and 71-^2 are exponential distributions 
with parameter 1), the high-SNR power offset is given according to Theorem [H by 

£ 00 = -max(E(log 2 |a| 2 ),E(log 2 |6| 2 )) = — - / e~ x log x dx = --L- (39) 

log 2 Jo log 2 

where the last equality is due to [28, pp. 567, formula 4.331.1]. Obviously this result coincides 

with the high-SNR power-offset derived by applying the definition of (see (fT3~T )) directly to 

the exact expression derived in [14] (see expression (1281)). 

Note that the same result holds if an attenuation factor is added to one of the fading paths, 

e.g. b m = ab m where b rn ~ £/V(0, 1) and a G [0, 1]; this follows directly from Theorem [Q but 

not from [14], which requires symmetric fading paths (i.e. a = 1). 

d) Intra-Cell TDMA and General Fading Statistic: Consider the following single user 
single-input single-output (SISO) flat fading channel for an arbitrary time index 

y = ax + n , (40) 

where x is the input signal x ~ CjV(0, P), and n is the additive circularly symmetric Gaussian 
noise n ~ CAf(0, 1). In addition, a is the fading coefficient a ~ n a satisfying conditions 
(HI). . . (H3) and known only to the receiver (receiver CSI). Assuming that the fading process is 
also ergodic in the time domain, the ergodic capacity of the channel is given by [24] 

C = E^log(l + P|a| 2 ) , (41) 

where the expectation is taken over the fading distribution ir a . Accordingly, under the mild 
conditions (HI). . .(H3), the high-SNR regime of this channel is characterized by 

<Soo = 1 ; Ax, = log 2 \a\ 2 . (42) 

Using Theorem [H we can now establish the following analogy between the multi-cell setup and 
the SISO channel at hand. 
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Corollary 4 The high-SNR characterization of the intra-cell TDMA per-cell sum-rate supported 
by the "soft-handoff" setup with fading distributions TT a , tci, such that K Wa log 2 \a\ 2 > log 2 \b\ 2 , 
coincides with those of a scalar single-user fading channel with fading distribution n a . 

This observation allows us to use the vast body of work done for the celebrated scalar flat 
fading channel [24]. In particular, the high-SNR characterization of flat fading channels with 
the following fading statistics have been considered in previous works: (a) Rayleigh distribution, 
(b) Rice distribution, (c) log-normal distribution, and (d) Nakagami distribution (see [24] and 
references therein). 

e) Intra-Cell TDMA and Opportunistic Scheduling: Throughout this work we have assumed 
that the instantaneous channel state information is known to the MCP receiver only. Here we 
further assume that some sort of ideal feedback channel is available between the MCP receiver 
and the K mobile users included in each cell. This feedback channel is used to schedule the 
"best" local user in each cell for transmission during the current time sloj^. In other words, in 
each cell the user with the strongest channel fade towards the BS located on the right boundary 
of each cell is scheduled for transmission^ with power P. Hence, the index of the selected user 
in the mth cell reads 

k m = argmax \a m)k \ 2 m = 1, 2, . . . , M . (43) 

k=l,2,...K 

The resulting M x (M + 1) channel transfer matrix Hm of this scheduling scheme is a two 
diagonal matrix with independent entries. The probability density function of the main diagonal 
i.i.d. entries' amplitudes is given by 

^,|a| 2 = K-k k ~ 2 X d^ , (44) 

following the maximum order statistics [30]. On the other hand, the i.i.d. entries of the second 
non-zero diagonal are distributed according to the original fading statistics n b . 

Assuming that ir K .12 and tx^ satisfy conditions (HI). . . (H3), we can apply Theorem [T|in order 
to derive the high-SNR characteristics of the per-cell sum-rate achievable by this opportunistic 
scheduling 

<Soo = 1 ; A*. = - max \^n KMa (log 2 y) , E^i, (log 2 |6| 2 )) . (45) 

4 See [29] for a similar scheduling deployed in the Wyner cellular uplink channel. 

5 Since the right most cell indexed (M+l), has no BS on its right boundary it randomly schedules a user for transmission. 
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For Rayleigh fading channels and in the case where the number of users per-cell is large K ^> 1, 
we can use the well known fact that the square of the maximum of the K amplitudes behaves like 
\ogK with high-probability (see [31]). Hence, the rate high-SNR power offset of this scheme is 

A^-logslogK , (46) 

revealing a multi-user diversity gain of log log K. It is noted that allowing additional power 
control to this scheme will yield better performances. However, we are unable to apply Theorem 
CD for this situation. Finally, choosing the BS located on the right boundary of the cell is arbitrary; 
taken the BS located on the left boundary of the cell yields the same results. 

V. Concluding Remarks 

In this paper we study the high-SNR characterization of the per-cell sum-rate capacity of 
the "soft-handoff ' uplink cellular channel with multi-cell processing. Taking advantage of the 
special topology induced by the setup, the problem reduces to the study of the spectrum of 
certain large random Hermitian Jacobi matrices. For the intra-cell TDMA protocol where only 
one user is active simultaneously per-cell we provide an exact closed form expression for the 
per-cell sum-rate high-SNR power offset for rather general fading distribution. Examining the 
result, it is concluded that in the high-SNR regime, the rate of the cellular setup at hand is 
equivalent to the one of a single user SISO channel with similar fading statistics. 

Turning to the capacity achieving WB protocol, where all K users are active simultaneously 
in each cell, we derive a series of lower and upper bounds to the rate. These bounds are shown 
(via Monte-Carlo simulations) to be tighter than previously known bounds. 

Note that in Theorem [2] points a) and c) and in Proposition |3j we take the fading coefficients 
relative to the users of one cell to be independent. Those results continue to be true if we 
assume correlation between the fading coefficients relative to the users of the same cell (but 
independence between cells). The proof is identical to the proof given in the paper. 

Some of the analysis reported here can be extended to include the case where Gm is (2p — 1)- 
diagonal for some p > 2 (e.g. p = 3 for the channel matrix of the Wyner model), using an 
adaptation of the "Thouless formula for the strip" derived originally in [32]. Using this approach, 
bounds similar to those of Prop. [3] may be provided on the rate. Details will appear elsewhere 
[33]. 
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Appendix 

A. Proof of Theorem [7] 

In order to streamline the proof we somewhat modify notation. We consider two random 
sequences of complex numbers (a n ) and (b n ). The (a n ) (resp. (b n )) are i.i.d of law 7r a (resp. 
lib) and the (a n ) are independent of the (b n ). We set Q = ((a n ), (b n )). We denote by P the 
probability associated with those random sequences and by E the associated expectation. For a 
given integer n, we consider a channel transfer matrix Hm of size M x (M + 1). 



H 



M 



ax 6i 
"•• 







a M 






b M 



\ 



We consider the following variable 

C M (P) = tr {log (/ + PH M Hl) } 



Note that, 



H mH 



a 2 b[ 




a\bi 

I |2 , |. |2 

|o 2 | + |o 2 | 







alb 





M^M-l 



\ 



t 1 2 1 2 

ffiMV-l \ a M\ + Pm| j 
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With this notation, as explained in Section III-Dl Theorem Q] follows from the following. 

Theorem 5 [K = 1] Assume (HI) and (H2). 

a) For every p > 0, Cm{P) converges F-a.s as M goes to infinity. We call the limit C(P). 

b) Further assume [(H3) or (H3')]. As p goes to infinity, 

C(P) = logP + 2 max (E Wa log |x| ; log \x\) + o(l). 

Proof of Theorem \5\ Without loss of generality, in the proof we can assume 
(H5) E Wa log|x| <E^log|x|. 

Indeed, we may exchange the role of entries a, and b{ for 1 < i < M by a right-left reflection, 
namely the transformation a,- = bM-j+i, % = a M -j+i, 1 < j < M. 

For part a), only (HI) and (H2) are needed. Since part a) is a consequence of general facts 
concerning products of random matrices and does not use much of the special structure in the 
problem, we bring it in Appendix ITJl 

Part b) uses the theory of Markov chains and is specific to the particular matrix H M . We note 
that as a by product of this approach, we obtain a second proof of part a), however under the 
additional assumption [(H3) or (H3')]. We provide a proof of Theorem [5] under the assumptions 
(HI), (H2) and [(H3) or (H3')] in Appendices E] and E 

The structure of the proof is as follows. We first introduce an auxiliary sequence which allows 
us to reformulate the problem in terms of a special Markov chain. The study of the latter, which 
forms the bulk of the proof of Theorem [51 is carried out in Section |B] 

1 ) Auxiliary sequence: We begin with a technical lemma. 

Lemma 6 Assume (H2). F-a.s, HmH\ 4 does not have multiple eigenvalues. 

Proof: We let D denote the discriminant of HmH\ { , it is a polynomial in 
{|aj| 2 + \bi\ 2 ,a i+ ib\,al +l bi} which vanishes when there is a multiple eigenvalue. Therefore, it 
is a polynomial in 9?dj, Saj, and ^bi It is not identically because for hi = and cij = i, 
the eigenvalues of HmH\ { are distinct. The result follows directly from the following lemma 
which is an easy consequence of Fubini's theorem. ■ 
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Lemma 7 Let Q be a function from C n to C We assume that Q is not identically and that 
Q(z\, . . . , z n ) is a polynomial in the and the Szj. Then the set of the roots of Q has Lebesgue 
measure 0. 

In the sequel, we denote by Ai > . . . > Xm the ordered eigenvalues of HmH' m . For a given 
A, we consider the following sequence (indexed by n) of complex numbers (the dependence in 
A will only be mentioned when it is relevant): x = 0, X\ = 1, and for n > 1, 

that is 

_ A — |a ra | 2 — |fr ra | 2 ^ Qn^j-i , 

Note that x A / + i(A) = if and only if A is an eigenvalue of H M H\ 4 . Moreover, x n+1 is a 
polynomial in A of degree n with highest coefficient 1/ nr=i( a l+:A)- One can tnus wr ite using 
Lemma [6] 

n n 

^ + i(A) = n^f+Ar 1 ri( A - p - a - s ' 

i=l i=l 

Hence, for A = — 1/P, 

1 1 M 

C M (P) = log(P) + — log |x M+1 (A)| + — lo S l^+i^l P " a - s - (48) 

i=i 

By the Law of Large Numbers (LLN), 

M 



lim "77 l0 S \ a i+l b i\ = lo g M + E vr 6 log |x| P - a.S. 

M— >oo ik/ ' * 

Because of (|48l , to prove Theorem [51 we only need to show the following lemma. 



Lemma 8 Assume (HI), (H2) and [(H3) or (H3')] 

a) For every A < 0, - log |x n+ i(A) | converges ~P-a.s as n goes to infinity. The limit is 7(A), 
the Lyapunov exponent defined by d62l ). 

b) Assume further (H5). Then 7(A) converges to E^ log |x| — E^ log \x\ as A goes to 0. 
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2) Reduction to a Markov chain: To prove Lemma[H we take c n = x n /x n -\, for n > 2. Note 
that by (l47l) and (H2), P-a.s, x n ^ 0, hence c n is well defined and non-zero. By (|47T) . we get 

A — \a n \ — \b n \ a n b n _i 



c n +i 



Let d n = c n alb n -i. Then, 



i \ I 1 2 I j, 1 2 l fl n| |&n-l| , ,2 i |2 / -, . |&n-l| 

d n+ i = A-|a n | - |6 n | = A - \b n \ - \a n \ | 1 + 



4 



Let e n = ( 1 + ^fJ- ) . Then d n+x = \ -\b n \ z - \a n \ z e n , and 



with the initial conditions, 



-A + 


a n-l 


2 


1 


-A + |6 n _i 


2 +l 


a n-l 


2 

e n -i 



(49) 



C-2 



A - |ai| 2 - |i>i| 2 



4&i 



c? 2 = A — |6i| 2 — | a-i | 2 . 

cf 2 G R and d 2 < — | £»i | 2 , hence, < e 2 < 1. From (|49b we conclude that for all n, e n e R and 
< e n < 1. Now, for all n, 



d n 1 



Then, 



-log X„ + i 
n n 



-, n+l 



8=2 
n+l 



(50) 



n 



i=2 



Hog 


6i-i 







- log(l - Ci 



converges to E,^ log \x\ — E Ta log |x| by the LLN. We now study in details the 



Markov chain e„. 



B. Study of the Markov chain e n and proof of Lemma \S\ 

For simplicity, we write 5 = — A and we re-index the chain so that it starts from e . As in 

A -U I n 1 I ^ 

(51) 



5 + 




i2 


1 


5 + |6 n _i 


2 +l 


fln-1 


2 

e„-i 
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We denote by P eo the law of the sequence starting from e and by E eo the associated expectation. 

Proposition 9 Assume (H2) and [(H3) or (H3')]. The Markov chain e n has a unique stationary 
probability, say, /is and for s G L 1 ^), for every starting point eo G [0, 1], F eo -a.s, 



Proof: We start with two lemmas that will be proved later on. 

Lemma 10 For a, (3,5 G M + , we define the function (j> a ^ (we suppress 5 from the notation) 
such that for e G [0, 1] 



For any given e G [0, 1], we define the sequence (9 n (e)) by 9 = e and for n > 1, 9 n (e) = 
^a,/3(^n-i( e ))- Then, <fi a> p has exactly one fixed point in [0, 1], say K a ^, and 9 n (e) converges to 
K a,/3- Moreover, the convergence is uniform in the starting point in the following sense: 



Finally if ol\ < (Xi and (3\ > (32, then K aij p ± < K a2t p 2 . 

Lemma 11 Assume (H2) and [(H3) or (H3')]. 

a) For eo G [0, 1], there exist two sequences (^(eo)) and (6^(eo)) in [0, 1] such that the law 
of e n under P eo and the Lebesgue-measure on [(^(e )), (^(eo))] are mutually absolutely 
continuous. 

b) (^(eo)) and (6^(eo)) converge to, say 6 1 and 2 respectively, 1 and 2 are independent 
of eo and Q 1 < @ 2 . Finally, the convergence is uniform in the starting point in the sense 
of Lemma UM 

c) If e G [G^G 2 ], then for all n, the law of e n under P eo is absolutely continuous with 
respect to the Lebesgue-measure on [Q 1 , G 2 ]. 

We recall some definitions from the theory of Harris Markov chains, which will be used 
extensively in the proof. We refer the reader to [34] for the relevant background. 




(Ve > 0)(3n G N)(Ve G [0, l])(Vn > n )(|^ n (e) - Ka , p \ < e). 
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Definition 12 Denote by (r n ) a Markov chain on I an interval o/R. Set I a probability measure 
on I, it is an irreducibility measure if for all measurable set A such that 1(A) > and for all 

r el 



l is a maximal irreducibility measure if it satisfies the following conditions: 

• I is an irreducibility measure. 

• For any other irreducibility measure V, V is absolutely continuous with respect to I. 



Definition 13 Denote by (r n ) a Markov chain on I an interval of R. A set A is called Harris 
recurrent if for all tq e A, F ro -a.s, the chain r n visits A an infinite number of times. The chain 
(r n ) is called Harris recurrent if given a maximal irreducibility measure I, every measurable set 
A such that 1(A) > is Harris recurrent. 

Definition 14 Denote by (r n ) a Markov chain on I an interval of R Denote by I a maximal 
irreducibility measure. For every measurable set A such that 1(A) > we denote by ta the time 
when the chain (r n ) enters A. A measurable set B is called regular if for every measurable set 
A such that 1(A) > 0, 



Definition 15 Denote by (r n ) a Markov chain on I an interval of R. Denote by A and B two 
measurable sets. We say that B is uniformly accessible from A if there exists an e > such that 



We continue with the proof of Proposition [9l Denote by I the Lebesgue-measure on [9 1 , 2 ]. 
By [34, Theorem 17.0.1], it is enough to prove that the Markov chain e n is /-irreducible, positive 
Harris with invariant probability fis- Denote B + the set of Lebesgue-measurable subsets of [0, 1] 
with positive /-measure. Here is a technical lemma that will be proved later on. 



(3n) P ro (r„6i)>0. 



. If 1(A) = then l{r : (3n) P ro (r„ G A) > 0} = 0. 
• For any irreducibility measure I', I is equivalent to 




sup E ro (r^) < oo. 



inf P ro ((3n) r n 6 



B) > e. 
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Lemma 16 Assume (H2) and [(H3) or (H3')]. For all B G B + , there exists n Q = n (B) such 
that for all n > hq, 

p n 4 inf P eo (e n G B) > 0. 

e 6[0,l] 

We continue with the proof of Proposition [9j 

Step 1: The Markov chain e n is l-irreducible, Harris and admits an invariant measure unique 
up to a constant multiple. By Lemma [T6l for eo G [0, 1] and B G B + , the chain has a positive 
probability to reach B in n steps starting from e . Therefore, the Markov chain e n is /-irreducible 
and by Lemma [EUc), I is a maximal irreducibility measure for the chain e n . For a given B E B + , 
by Lemma [HI the chain e n has a probability at least p no to reach B in n steps, hence the chain 
will eventually reach B and hence come back to B an infinite number of times, therefore B is 
Harris-recurrent and the Markov chain e n is Harris. By [34, Theorem 10.0.1], the Markov chain 
e n admits an invariant measure unique up to a constant multiple. 

Step 2: The Markov chain e n is aperiodic. By [34, Theorem 5.4.4], there exists an integer d, 
the period of the chain, such that there exist disjoint measurable sets D , . . . , Dd-i such that 
• For i = ... d — 1, if 6i E Di, then P ei (ej + i G A+i) = 1 (mod d). 

. I ((uf =1 A) c ) = 0. 

By Lemma [TH for ni > n large enough and n > n\, the Lebesgue-measure on J = [(26 1 + 
9 2 )/3, (O 1 + 20 2 )/3] is absolutely continuous with respect to the law of e n under P eo . Therefore, 
for any n > m, if e n G Di, then J C -Dj, and then, if d > 1, e n+ i G D i+ % and thus also J C A+i, 
a contradiction. Hence, d = 1. 

Step 3: The set [0, 1] is regular for the Markov chain e n . Take B G B + . By Lemma [T6l the 
time it will take for the chain e n to enter B is a.s bounded above by n times a geometric 
random variable of parameter p no , hence it expectation is bounded above by n /p no , hence [0, 1] 
is regular. 

Now we apply [34, Theorem 13.0.1] and get that the Markov chain e n is positive Harris, 
hence has a unique invariant probability that we denote fis- 

U 

Proof: [Proof of Lemma [T6ll 
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The Lebesgue-measure on [O 1 , 2 ] is regular hence there exists an e > such that B fl [0 1 + 
e, G 2 — e] has positive Lebesgue-measure. By Lemma [TT] a) and b), we can take n such that 
for any given n > n and any given starting point e , P eo (e n G B) > 0. Fix n > no. Set 
■0(e o ) = P eo ( e n & B). By (H2), ^ is a continuous function on [0, 1]. By compactness, 

inf P eo (e n G B) > 0. 
e e[o,i] 

■ 

Proof: [Proof of Lemma [TTIl Let us start assuming (H3'). 

a) We first assume that Ai a ,Aib G We use the notation of Lemma ITOl For eo G [0, 1] 
and n, we define 0*(e ) = 0m a ,/v! b ( e o) an d 0^(e ) = 0vU a m(l ( e o)> where n is the n-th iteration 
of the function <r>. Note that for el < e 2 G [0, 1], al < a 2 G K + and /51 < /3 2 G K + , 

V> : [ei,e 2 ] x [a 1; a 2 ] x [(3 U (3 2 ] — ► [0 Q1 ,/3 2 (ei), a2 , A (e 2 )] 

(x,a,(3) i — > 0a,/3(e) 

is well defined and onto and the inverse image of an interval which is not a singleton has positive 
Lebesgue-measure. Therefore, by induction, the Lebesgue-measure on [0*, 6* 2 ] is absolutely 
continuous with respect to the law of e n under P eo . Moreover, by (H2) and (I5TI) . the Lebesgue- 
measure on [0*, &n) an ^ me ^ aw °f e n under P eo are mutually absolutely continuous. 

b) It is a direct consequence of Lemma ITOl and we get 1 = K majMb and 2 = KM a ,m b - By 
Lemma [TO] and (H3'), n maMb < K MatTnb , hence 1 < 6 2 . 

c) 4>m a ,Mb * s increasing and K ma:Mb a fixed point hence if K ma) M b < e , then for all n, 
«m a ,A^ 6 < ^i(eo)- In the same way, for all n, KM a ,m t > C(e )- 

If .M a = oo (resp. Mb = oo), we take for all n > 1, 6> 2 = 1 (resp. 0* = 0) and 6 2 = 1 (resp. 
1 = 0) and the proof is the same. 

Let us now assume (H3). The proof is the same with for all n > 1 and all eo G [0, 1], 
0i(e o ) = 0, for all 72 > 1 and all e G [0, 1] (except for n = 1 and e = 0), # 2 (e ) = 0. We get 
1 = and 2 = 1. ■ 
Proof: [Proof of Lemma EJ For e G [0, 1], 

- {6 + (3 + ae y 

<j)' a g is decreasing and <f>' a Jl) < 1. If </>' a *(0) < 1, then Q /3 is contracting hence admits a fixed 
point and its iteration on any starting point converges to the fixed point. Suppose 0^(0) > I. 
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Denote by e the only point of [0, 1] such that (f>' a n(e) = 1. Set <p a ,p(e) = 4>{e) a ,p — e. Then 
0a,/?(O) > 0, 0a 1( a(l) < 0, and Qj/ g is increasing on [0, e] and decreasing on [e, 1]. Hence, 
0a,/3(e) > and <fi aj p is on exactly one point which is a fixed point for 4>a,p- We denote that 
fixed point K ai p. If e G [K, a< p, 1], since 0^ is increasing, for all n, 9 n (e) G [K a ,p> 1] and aiJ g is 
contracting on [k Qj/3 , 1] hence 0„(e) converges to If e G [0, K a} p], for all n, 9 n (e) G [0, K a ^], 
and (j) a ^ is non-negative on that interval, hence 9 n (e) is non-decreasing. Therefore, it converges 
and since a /3 is continuous, the only possible limit is K a> p. To prove the uniformity in the 
starting point, we use the fact that <p a ^ is increasing, hence for all e G [0, 1] and n, 

9 n (0) < 9 n (e) < n (l). 

That gives the uniformity. Finally, assume a\ < ot2 and Pi > fa- 4>a,p(e) is non-decreasing in 
a, decreasing in (3 and non-decreasing in e hence by induction, 0^ ^(0) < 0™ 2 p 2 (0), where <\> n 
is the n-th iteration of the function 0. Hence, K a i,Pi < K a 2 ,p 2 - If K ai,/3i = K a 2 ,p 2 , men 

K ai,/3i — 0ai,/3i ( K oi,/3i) < 4>a 2 ,p 2 ( K oi ,/3i ) = 4>a 2 ,f3 2 ( K c*2,/32 ) = K a 2 ,p 2 i 

which gives a contradiction. 

■ 

We continue with the proof of Lemma [8l Recall that < e n < 1, hence fis is stochastically 
dominated by an atom at 1. fig is the invariant measure, since the function (f) a ,p{-) is increasing 
in e, fi$ is stochastically dominated by the law of the chain started at 1 after one step: 



■< £ , . „ ,2 . , ,2 =< C 



5 + 


\ao\ 


2 


5+\ 


bo\ 


2 + 




2 





|«o| 


2 




2 + 




2 



I |2 

Thus, denoting by tt the law of 2 , and using (HI), 



y — log(l — x)dfj,s(x) < J — log(l — x)rf7T ( 



X < oo. 



That is 

-log(l--) GL 1 ^). (52) 

With Proposition HI we get 

- V]-log(l-eji.) >/ -log(l - x)d/i s (x) P £2 -a.s. (53) 

k=2 J{i 
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With (|5Q|) . it gives a proof of Lemma [8] a). 

Let us prove Lemma [8] b). Take r] > and e > small. 

i 

-log(l - x)dfi S {x) 







= / -log(l - x)dfi 5 (x) + / -log(l -x)d^ s (x) + / -log(l -x)d^ 5 (x) (54) 

< -elog(l-e) -bg^i([e, 1]) + / - log(l - x)dp, s {x). 

Ji-n 

By (|52l) . the last term converges to as 77 goes to 0. By (l50l) . (l53l) and (l54l) . to prove Lemma [8] 
b), we only have to prove that for any given e > 0, 

1])— >0. 

0— HJ 

For that, by Proposition HI we need to show that the proportion of the time that the chain e n 
spends above e converges to as 5 goes to 0. We take < e < e < 1, where e will be 
chosen later. We consider the Markov chain z n = loge n and the random function g n such that 
z n = 9n{z n -i). It is enough to show that the proportion of the time that z n spends above \oge 
goes to as 5 goes to 0. Let us couple z n with another Markov chain w n , such that w n > z n 
a.s. and that the proportion of the time that w n spends above loge goes to as S goes to 0. 
For that, we need good information on the jumps of z n . 

Lemma 17 Assume (HI) and (H5). Set 

jnK^n—x) Z n Z n — 1 

= log (^7 + \a n -i\ 2 ^ - log (5+ |6„_i| 2 + \a n ^\ 2 e^- 1 ) . 

(V5 > 0) (Be' > 0) (Vx > logs') 

a) Ej n (x) < 0, 

b) Var Jn (x) <V = E ((logda^^ 2 + l^ 2 )) 2 + (logda^i) 2 )) 2 ) + C. 

C is a constant independent of everything, e' is a function of 5 but we will not write it to keep 
the notation clear. Moreover, 

lime' = 0. 

The proof will be done at the end of the section. 
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We continue with the proof of Lemma [8] b). We take 5 > such that < e' < e < e$ < 1. 
We define w n in a way that it stays between logs' and 0. Set u> = z , for 5 small enough, 
w > \oge'. For x G [logs'; 0], denote 



h n (x) = g n (x) - Ej n (x) > g n (x). 



That is 



h n (x) = x + lot 



- + 




-1 


2 


5+\ 


K-i\ 


2 + 




2 e x 



Elog 



- + 




-i 


2 


^ + l&n-ll 


2 


+ 


a n-l\ 


2 e x 



(55) 



Note that 



E(h n (z n -i) - Zn-l|*r»-l) = 0. (56) 



• If h n (w n -i) > 0, set w„ = 0. 

. If ^(iDn-i) < logs', set u> n = logs'. 

• Otherwise, set w n = h n {w n -\). 

In the first two case, we say that the chain is truncated. Note that for all n, w n > z n . Indeed, 
either w n = > z n or w n > h n (w n -i) > g n {w n ^) > g n {z n -i) = z n , by induction and using 
the fact that g n is a.s non-decreasing. Therefore, the proportion of the time that the chain w n 
spends above logs is larger that the proportion of the time that chain z n spends above logs. 

Proposition 18 Assume (H2). 

a) The Markov chain w n has a unique stationary probability, say, v§ and for s G L 1 ^), for 
every starting point wq G [logs', 0], F WQ -a.s, 



-Y^ s(wi) > / sdus. 

i=0 J 



b) We denote T the return time to 0, starting from 0. Then ^(0) = 1/E T. 



Proof: See [34] and Definitions [1211151 for the theory of Harris Markov chains that we will 
use extensively in the proof. Define the following probability measure on [logs', 0]. For B a 
Borel set, 

n=0 
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Let us prove that the Markov chain w n is I-irreducible, positive Harris with invariant probability 
vs. By [34, Theorem 17.0.1], that will prove a). We use the following lemma that will be proved 
later on. 

Lemma 19 Assume (H2). 

a) There exist c > and 9 > such that for all x E [logs'; 0], 

POnO) >X + C)> 9. 



b) Set N 



log e' 



. is a recurrent point for the chain w n and the time between two visits 
at is a.s bounded above by N times a geometric random variable of parameter 9 N . 



We continue with the proof of Lemma [181 The sets which have positive /-measure are exactly 
the sets that have a positive probability to be visited starting from 0. Moreover is a recurrent 
point. Therefore, the Markov chain w n is I-irreducible and 7 is a maximal irreducibility measure. 
Moreover, take B with positive /-measure, B is uniformly accessible from {0}. Therefore, we 
can apply [34, Theorem 9.1.3 (i)] and since is Harris-recurrent, B is also Harris-recurrent, 
therefore, the chain w n is Harris-recurrent. By Lemma [19] b), the time between two visits at 
has finite expectation (bounded above by N/9 n ). Therefore, by [34, Theorem 10.2.2], the chain 
w n is positive-Harris and admits a unique invariant probability measure. That finishes the proof 
of point a). The point b) is a consequence of 

l = ^([logs',0]) = ^(0)E [T], 

which comes from [34, Theorem 10.0.1], which we apply to A = {0}, which has positive 
/-measure. ■ 
Proof: [Proof of Lemma [T9l a) We consider here 5 > 0. We denote by Supp(X) the support 
of the law of a random variable X. We take 5 small enough. We consider for x E [logs'; 0] the 
function 

<t>(x) = max{y; y G Supp(/i n (x) - x)}, 

which by (H2) and (|55l) is a continuous function of x. Moreover, since ~E(h n (x) — x) = 0, is 
strictly positive. By compactness, there exists c > such that for x E [logs'; 0], 

6(x) > 2c, 
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F(h n (x) >x + c)>0. 

By (H2) and (1551) . F(h n (x) > x + c) is continuous and once again, by compactness, there exists 
9 > such that for x G [loge:'; 0] , 

F(h n (x) >x + c)>6. 

b) If there are at least N steps in a row such that h n {w n -i) > x + c, then the chain reaches 0. 
By the point a), that happens with probability at least 9 N > 0, hence is a recurrent point for 
the chain w n and the time between two visits at is a.s bounded above by N times a geometric 
random variable of parameter 9 N . ■ 

We continue with the proof of Lemma [8] b). By Proposition [18] a), to prove that the proportion 
of the time that w n spends above \oge goes to as 8 goes to 0, we only need to prove that 

^([log £ ,0])— >0. 

Let us first prove that ET - — ► oo, which by Proposition [TBI b) will prove that 

u s (0) ► 0. 

We use the following lemma. 

Lemma 20 Assume (H2). 

a) There exist u > and a > dependent on e and independent of 8 such that for all 
x e [2\oge; 0], 

P (h n (x) > x + u) > a. 

b) There exist v > and (3 > dependent on e and independent of 8 such that 

F(\oge < hi(0) < -v) > (3. 

The lemma will be proved later on. 

We continue with the proof of Lemma [8] b). We denote A the event logs < h\(ti) < —v. On 
A, we define the stopping time 

f = 1 + inf{n > 1; h n+1 (w n ) > or h n+1 (w n ) < loge'}. 
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We now condition on the event A and on x = ^i(O), denote by P and E the associated 
probability and expectation. T < T is the first time the chain is truncated. Moreover, for n <T, 
Wn = h n (w n -i), so with (|56l) . by classical martingale arguments, 

E (/if (wf_i)) = x . 

We denote by A$ the event that w n reaches logs' before 0, we set p = F(Ao), X = 
E (hf (wf_i)|>lo) and Xi = E (hfiwf^lAo). 

x = pX x + (1 - p)X . 

_ Xq — Xq 

Xq — X\ 

X > and X\ < log e' < log e < x < — v hence, 

. -Xq V 

**=x x *=x x - (57) 

Using X\ < logs' , (f57l) and w 2 n — Vn, which is a super-martingale by Lemma [T71 b). 

E(T) > E(T) > _LJ: V T - i;y 



> 



> 



> 



> 



V 

pXl + (l-p)Xl-xl 
V 

pXj - xl 
V 

f(-Xi) -x\ 
V 

v(— logs') — Xq 

V ' 



We integrate over x and use F(A) > (3 and E(/ii(0) 2 |^4) < (logs) 2 . 

v{-loge') -E(/i!(0)V) 



E(T|^) > 



E(T) > /3 



r 

f ( — logs') — (logs) 2 



We have proved that ET ► oo, which proves that za(0) ► 0. 

<5^0 <5->0 

Using Lemma l20l and the invariance of v$, let us prove by induction that for M < |" ~ ° g£ ] , 

vs{[-Mu]0])<a-"v s {Q). 
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u s ([-(M - l)u; 0]) > J Us(dwo)W m ( Wl e [-(M - IK 0]) 

> / u 5 (dw )W W0 ( Wl £ [-(M - l)u; 0]) 

> / u s (dw )¥ Wo (hi(w ) > u + w ) 

> ais s ([-Afu;0}). 

Therefore, 

i/ 5 ([logs; 0]) < a\^^u s (0). 

So, 

„,(po ge; 0]) — >0. 

That concludes the proof of Lemma [8] b). 

Proof: [Proof of Lemma l20ll We consider here 5 > 0. We denote by Supp(X) the support 
of the law of a random variable X. We take 5o small enough, 
a) We consider for x E [2 loge; 0] and < 5 < S the function 

4>(x, 5) = max{y; y G Supp(/i n (x) — x)}, 

which by (H2) is a continuous function of (x, 5) because (h n (x) — x) is continuous in (x,S). 
Moreover, since M(h n (x) — x) = 0, is strictly positive. By compactness, there exists u > 
such that for x £ [2 \oge; 0] and < 5 < 5 , 

4>(x, 5) > 2u, 

F(hn(x) >x + u)>0. 

By (H2), ¥(h n (x) > x + u) is continuous and once again, by compactness, there exists a > 
such that for x £ [2 \oge; 0] and < 5 < 5 , 

¥{h n {x) > x + u) > a. 

b) For all < S < S , there exist e > and v > such that P(log£:o < ^i(0) < —v) > 0. Like 
in the proof of a, by (H2), we can chose e > and v > continuous in 5. By compactness, we 
can chose Sq > and t> > independent of 5 such that for all < 5 < So, P(log£o < ^i(0) < 
— v) > and like in the proof of a), by (H2), that probability can be chosen continuous in 5. 
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Therefore, by compactness again, there exists (3 > dependent on e and independent of 5 such 
that P(logs < hi(0) < -v) > (3. Take e < e , we have P(logs < hx(0) < -v) > (3. U 
Proof: [Proof of Lemma [T71 Note that by (HI), V < oo. j n (x) is a non-increasing continuous 
function of x and so is Ej n (x). Ej n (0) < 0, hence given 5, there exist < e' < 1 such that 
Ejn (logs') < 0, and for x > logs', Ej n (x) < 0. That gives point 1. For point 2, take C such 
that for all x > 0, 

(log(ac + l)) 2 < (log(x)) 2 + C. 

To prove that lim^^o s' = 0, it is enough to prove that for a given L < 0, we can find <5 small 
enough such that Ej n (L) < 0. That is true because for a given L, Ej n (x) is a continuous function 
of 8 which, by (H4) is negative for 5 = 0. ■ 



C. Proof of Theorem [2] 

We reformulate the problem in the spirit of Appendix lAl Let K > 1. The a, L (resp. 6j) are 
now independent complex vectors of size K whose coefficients are independent and distributed 
according to ix a (resp. 7Tb). We denote by P the probability associated with those random 
sequences and by E the associated expectation. We consider the following M x K(M + 1) 
channel transfer matrix: 



H 



M 



'■• 











\ 







V 



a M b 



We consider the following variable 



C M (P) = —tr{\ot 



P 
K~ 



where P = Kp. Note that, 

^1 |2 i |l |2 
| CZi | + |Oi| 

< a 2 ; &! >t 



HmH m 



< a 2 ; b x > 



I+-H M Hl 







\a 2 \ 



\b 2 \ <a 3 ;b 2 > 







< a M ] &M-1 > 





< a M', bu-i > 
I |2 . it |2 

% + \o M \ 
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where |aj| = J2k=i K&l and < a u bj >= J2k=i( a i,kVb jjk . 

Theorem 21 Assume (HI), (H2) and (H4) 

a) For every p > 0, Cm{P) converges F-a.s as M goes to infinity. We call the limit C(P). 

b) As P goes to infinity, 

C(P)=logP + Elog (^^~\ +o(l), 

where the expectation is taken in the following way. e and b are independent, b is a complex 
K-vector whose coefficients are independent and distributed according to tti,. The law of 
e is too, which is the unique invariant probability of the Markov chain defined by 

e n + |6 n -i| 2 sin 2 (a n ,b n _ i ; 



I I 2 



e„ + |&n-i| 2 



The rest of this appendix is devoted to the proof of Theorem | 
As in Appendix lAl we define the sequence x n as follows. x = 0, x\ = 1, and for n > 1, 

_ A - |a„| 2 - |b n | 2 < a n ; b n _i > f 

Xn+l — T X n T x n-l- \Po) 

< a n+l] On > < a n+l] °n > 



We get, like in <g8]), for A = -1/p, 

1 1 - 

C M (P) = \og(P/K) + — log |x M+ i(A)| + — lo S l< a ^ °i>\ P " a - s - ( 59 ) 



i=i 

Set c n = x n /x n -i, for n > 2. By (|58l) . we get 

A - \a n \ 2 - \b n \ 2 <a n ;6 n _i> t 



< a>n+i] b n > c n < a n+ i; b n > 
Let d n = c n < a n ; b n _i >. Then, 



d n +i — A — \a n \ — \b n 



I » |2 

2 it 1 2 l< a ni °n-l >\ 



, ,2 i ,2/, , l b n-i| cos 2 (a n , 6 n _i^ 
A - |6J - |aj I 1 + 



dn 



where 

cos 2 (a„, 6 n _ x ) = |< a n ; b n _i >| 2 / |a n | 2 |6„_i| 2 . 

Note that < cos 2 < 1. Let e n = —d n — |6 n _i| 2 . 

, , e n + |6 n _i| 2 sin 2 (a n , b n _i) . 
e„+i = -A + (a^l ( -J ) , (60) 

e n + On-l 
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where sin 2 = 1 — cos 2 . With the initial conditions, d 2 < — \bi\ 2 , hence e 2 > and for all n, 
e n > 0. Note that (e„) is a Markov chain and that for all n, e n is independent of a n and b n -i- 
By (EH), we get 

^ M+1 A/ 

C M (P) = log(P/if) + — ]T log |c,(A)| + — ^ log |< a i+1 ; 6, >| 

i=2 «=1 
1 M+1 

= Iog(P/lf) + - MKI) + °(!) ^ 61 ) 

= ,„ g( P ) + 1 Lg,o g (M^d!j +o(1) 

We only need to study the Markov chain (e n , 6 n _i). For convenience, we set 5 = —A and we 
allow 5 = 0. We also assume without loss of generality that the chain starts at (ei, &o)- 

Proposition 22 Assume (H2) and (H4). Take 5 > 0. The Markov chain (e n (5),6 n _i) /za* a 
unique stationary probability, say, fi$ and for s G /or every starting point (ei,b ) G 

R + x C^, P (ei ,i, o) -a.5, 




Moreover, p,g is weakly continuous in 6 = 0. 

Proof: We consider the Markov chain (e„) on the compact [0, oo]. By (l60l) . for n > 1 and 
e G [0, oo], P e (e n = oo) = 0. Consider (|60l ), by (H2), for ei G [0, oo), the law of e 2 under P ei is 
absolutely continuous with respect to the Lebesgue measure on [5, oo]. Moreover, by (H4), the 
law of 62 under P ei and the Lebesgue measure on [S, oo] are mutually absolutely continuous. 
Therefore, for e x G [0, oo) and n > 3, the law of e n under P ei and the Lebesgue measure 
on [S, oo] are mutually absolutely continuous. That fact allows us to prove like in Appendix 
|B] that the Markov chain (e„) is /-irreducible, positive Harris with invariant probability ms, 
where / is the Lebesgue measure on [5, oo]. Since P e (e n = oo) = 0, ms does not charge {oo}. 
We identify mg and the measure it induces on R + . We denote by n b the law of b. Since for 
n > 1, e n and b n _i are independent, the Markov chain (e n , fe n _i) is I x n b -irreducible, positive 
Harris with invariant probability fi$ = ms x lib. By [34, Theorem 17.0.1], the Markov chain 
(e n (A), has a unique stationary probability jis and for s G L 1 ^), for every starting point 
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( ei ,6o) GM+ x C K , P( ei ,b )-a.s, 

-Y^s(e i ,b < _i) ► / sd[i s - 

i=0 J 

Let us prove that ms converges weakly to mo when 5 converges to 0, which will finish the 
proof, {ms, 5 > 0} are measures on the compact [0, oo] hence it is enough to show that m 
is the only limit point when 5 goes to 0. By (H2), for a point x and an interval A in [0, oo], 
Pei(e2(<5) E A) converges to P ei (e2(0) E A). It implies that a limit point must be an invariant 
measure for the chain with 6 = 0. The only possibility is m . ■ 

By (l60l) . ms is stochastically dominated by the law of \a n \ 2 + 5. Therefore, by (HI), (x, y) — ► 
log(x + y) E L 1 ^)- (I6T1) and Proposition [22l conclude the proof of Theorem EU 

D. Product of random matrices 

We prove Lemma[8]a) assuming only (HI) and (H2). We use the theory of product of random 
matrices theory. For a general introduction to the aspects of the theory we use here, the reader 
may consult [25], [26], [35]-[37]. 

Let us take |-| any norm on C 2 and ||-|| the associated operator norm on matrices. For a given 

A-|an| 2 -|bn| 2 a " fe u-l 

1 

For a,a',b,b' E C — 0, we define the following invertible matrices 

A-|"| 2 -|fr'| 2 ofct 



/ Xn+1 ' 


H 







)i 











g(X, a, a', b, b 

Finally, we define 

( A-|an| 2 -|b„| 2 _anbl_ 1 
1 

M n = g n ...g v 

So that 



j x n+l ' 










\ ^ J 
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Set 8 = (C — O) 4 which is a Borel set of a separable and complete metric space. X n = 
(a n+ i, a n , 6 n , 6 n _i) is a Markov chain on £, with invariant measure II = 7r a x 7r a x ir b x 7r&. With 
(HI), 

E n (log+ \\g(X,a,a',b,b')\\ + log + ||g(A, a, a', 6, fr') - '!) < oo. 

Notice that gr n (A) is a continuous function of X n , therefore ((X n , M n ), II) is a multiplicative 
Markovian process. By [38, Example 1 and Proposition 2.5], 1/nlog ||M" n (A)|| converges P- 
almost surely and in L^fi), we set 

7 (A) = lim -log||M n (A)||. (62) 

n— >oo U 

7(A) is the first Lyapunov exponent. 

The L^fi) convergence already gives an easy upper bound for 7(A). By the property of 
operator norm, 

7 (A) <E n log IMA)!). 

Moreover, we can refine that bound into a whole family of upper bounds, for fceN, 

7(A) <^E n log ||^(A)...^(A)||. (63) 

Note that this upper bound is getting better as k increases and tight as k — > 00. 
Let us now prove that 

-log|x n+ i(A)| ► 7(A). 

fl n— +00 

Definition 23 The multiplicative system ((X n , M n ), U) is irreducible if there is no measurable 
non-random family {V(X),X G E} of proper subspaces of C 2 such that 

M n V(X ) = V(X n ), F-a.s, Vn G N. 

Lemma 24 Assume (H2). The multiplicative system ((X n , M n ), IT) is irreducible 

The proof is an adaptation of the proof of [39, Proposition 6.1.1]. 

Proof: The proof is by contradiction. Assume that there is a measurable family {V(X), X G 
E} of proper subspaces of C 2 such that 

g 3 V(X 2 ) = V(X 3 ), P-a.s., Vn G N. 
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We parameterize the proper subspaces of C 2 by ( ] for c in (—00, 00]. There is a measurable 



family {c(X),X G E} such that g 3 
computation gives 

0(04,03,63,62) : 

that is 

c(o 3 , a 2 , b 2 , W) 







' c(X 3 )\ 


| c(X 2 )) 


and | 











\ I 1 2 17 i2 

A — |a 3 | — J 03 1 
alb, 



a 3 b\ 



4^3 



c(a 3 , a 2 , 6 2 , 6i)4 fe 3 



are P-a.s. collinear. A direct 



P-a.s. 



a 3 b\ 



a\b 3 ( A |03 { |b31 0(04,03,63,62) 



P-a.s.. 



Note that the RHS does not depend on a 2 and 61, hence, c(a, a', 6, 6') does not depend on a' and 
6'. Setting rf(a, 6) = o)b c(a, 6), we get 



ci(a4, 63) = A — |a 3 



1 1 2 1 j, 1 2 

|2 it |2 a 3 °2 

! - N 



P-a.s.. 



d(a 3 ,6 2 ) 

The RHS does not depend on 04, hence, d(a, 6) does not depend on a. From (|64|) . we get 



(64) 



d(b 2 



03 



d(6 3 ) - A + |a 3 | 2 + |6 3 | 2 ' 



P T -a.s.. 



The RHS does not depend on 6 2 , hence, d(b)/ |6| 2 does not depend on 6, set d(b) = L \b\ 2 , where 



L is a fixed constant. Then, 



(L + l)|6 3 r = A-|a 3 r 1 + - , P.-a.s.. 



1 



1 1 2 

If L ^ —1, |6 3 | is a measurable function of a 3 and since it is also independent of a 3 , it is a 
constant, which is in contradiction with (H2). Hence L — — 1, which gives a contradiction with 
A < 0. ■ 
By [38, Lemma 2.6], irreducibility implies that 



lim — loj 



^n+2 



7 P-a.s. 



The following lemma completes the proof. 
Lemma 25 Assume (HI). 



lim — I loj 



Xn+2 



— log|x n+ i| =0 P — a.s. 



June 16, 2008 



DRAFT 



41 



Proof: 



log 




log \x n 



+ll 



log 




> 0. 



Let us prove that for e > 0, P I - log 




> e \ is a summable series, which by the Borel- 



Cantelli Lemma will prove the lemma. We have 



P I - log 

n 




> : I <P| ^log(|c n | + l) >e 



<P(M >e" £ -l) 
< PflcJ > e~) 



\b n -l 


1 




a n 


1 - e n 



< p 



< P ( J^zli > e ? 



> e 2 
fP 



> e 4 



(65) 



| CJj-^ | J \ 1 ^77, 

We analyze the right side of (l65l) . We use the fact that log \a n \ and log |& n _i| have a second 
moment by (HI) and that it does not depend onn. By the Bienayme-Tchebicheff inequality, we 
get 

ne s 
T 

J-OIC I t U)U i7„- I - |U|> |f/„ I 



P r 1 > e 4 



P(tog|6„_i| -log|a n | > 
16E((log|6 n _i| -log|a„|) 2 ) 



(66) 



n 2 e 2 



implying that the first term in the right side of (l65l) forms a summable series. Moreover 

1 



log- 



< log 



-A + |6 n _i| 2 + |a n _i| 2 



-l e n \b n -i\ 
which has a second moment by (HI), hence, by a computation like (|66l) and the Bienayme- 
Tchebicheff inequality, P (tz^t > e ^ J is a summable series. The Borel-Cantelli Lemma applied 
to the right side of (1631) concludes the proof. ■ 



£. Determinants of Jacobi Matrices 

An interesting and useful characterization of an M x M Jacobi matrix is that its determinant 
can be expressed by the following recursive formula [40] 

det G m = fxi,m dct G m —i [G m ] mim _i[(r m ] m _i jm det Gm—2 



; m = 3, . . . , M , (67) 
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with 

detGx = [G m ] 1A 

(68) 

det G2 = [G m ]i t i[G m ]2 t 2 — [G m ]i j 2[G ! m ]2 i i , 
where G m is the principle submatrix of Gm, obtained by deleting its last (M — m) columns. 
This characterization already used by Narula [14], can be easily proved by straight forward 
calculation of the determinant of Gm, starting from its last row. 

Examining (1671) . it is observed that the determinant of a square Jacobi matrix is dependent on 
a weighted sum of its two largest principle matrices' determinants only. Furthermore, det G m _i 
and detG m _ 2 are independent of the entries [G m ] OT , m , and [G m ] m+hm . 

It is worth mentioning that this approach can not be extended for matrices with a number of 
non-zero diagonal higher than 3. Hence, a similar formula, can not be obtained even for five- 
diagonal matrices and the resulting formula involves O(M) determinants of submatrices (not 
necessarily principle submatrices). 
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High-SNR Power Offset (K=2) 
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Fig. 2. High-SNR power offset bounds for Rayleigh fading, K — 2, and bounds order n = 1, 2 • • • ,8 
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High-SNR Power Offset (K=10) 
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Fig. 3. High-SNR power offset bounds for Rayleigh fading, K = 10, and bounds order n = 1, 2 ■ ■ ■ ,8 



High-SNR Power Offset (n=2) 



V — Upper bound 
A — Lower bound 

Previously known lower bound 

Asymptotic (K==°) 




5 6 7 

Number of Users K 



Fig. 4. High-SNR power offset bounds (order n = 2) for Rayleigh, and K = 2,3,4,6,8,10. Note that for K = 1, 
~ !og"2 ~ °- 833 - 
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